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SUMMARY 


A theoretical investigation has been made to determine the effect of 
nonlinear stability derivatives on the laterai stability of an airplane. 
Motions were calculated on the assumption that the directional- stability 
and the damping-in-yawing derivatives are functions of the angle of side- 
slip. The application of the Laplace transform to the calculation of an 
airplane motion when certain types of nonlinear derivatives are present 
is described in detail. The types of nonlinearities assumed correspond 
to the condition in which the values of the directional-stability and 
damping-in-yawing derivatives are zero for small angles of sideslip. 

The results of the investigation indicated that under certain condi- 
tions the nonlinear stability derivatives assumed in the analysis caused 
a motion which had different rates of damping for the large and small 
amplitudes of motion, with very little damping at the small amplitudes. 

In general, the period of the resultant oscillation increased with time. 


INTRODUCTION 


Recent flight tests of several airplanes designed for high-speed 
high-altitude flight have indicated neutrally damped lateral oscillations 
of small amplitude generally referred to as snaking. Upon examination of 
the flight records, the decrement of the oscillatory motion is found in 
some cases to be different for the large and small amplitudes of motion 
with a neutrally stable oscillation occurring at the small amplitudes. 

One of the explanations offered for the cause of this type of motion is 
that some of the stability derivatives are nonlinear ; that is, the 
derivatives have different values for the large and small amplitudes of 
motion. The nonlinearity could be caused by boundary-layer effects or 
flow separation due to poor fairing at the junction of the tail surfaces. 
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The present paper represents a preliminary investigation of the effect 
of the presence of two nonlinear stability derivatives, the directional- 
stability derivative C n p and the damping- in -yawing derivative C nr * 

on the motion of an airplane. These derivatives were selected for the 
analysis since the damping of the oscillation is a function of C nr and 

since C nr depends upon the C n p contributed by the tail. The 


derivatives C n p and were both assumed to be functions of the 

sideslip angle p. Calculations were made of the airplane motion due to 
a disturbance in sideslip for three different types of nonlinearities. 
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SYMBOLS AND C OEFFIC 3EM?S 

angle of roll, radians 
angle of yaw, radians 

angle of sideslip, radians except where noted in figures (v/v) 

yawing angular velocity, radians per second ( dty/dt ) 

rolling angular velocity, radians per second (d0/dt) 

sideslip velocity along the Y-axis, feet per second 

airspeed, feet per second 

mass density of air, slugs per cubic foot 

dynamic pressure, pounds per square foot 

wing span, feet 
wing area, square feet 
weight of airplane, pounds 
mass of airplane, slugs (W/g) 

acceleration due to gravity, feet per second per second 
relative density factor (m/pSb) 
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Sc 

" n p = dT 


n 


inclination of principal longitudinal axis of airplane with 
respect to flight path, positive when principal axis is 
above flight path at the nose, degrees 

angle of flight path to horizontal axis, positive in climb, 
degrees 

radius of gyration in roll about principal longitudinal 
axis, feet 

radius of gyration in yaw about principal vertical axis, 
feet 


nondimensional radius of gyration in roll about principal 
longitudinal axis (k Xo /b) 

nondimensional radius of gyration in yaw about principal 
vertical axis (kz 0 A*) 

nondimensional radius of gyration in roll about longitudinal 
stability axis (\/k Xo 2 cos 2 t) + Kz 0 2 sin 2 T]^ 


nondimensional radius of gyration in yaw about vertical 
stability axis ^Kz q 2 cos 2 t] + K Xo 2 sin 2 T]j 
nondimensional product-of- inertia parameter 
^ K Z 0 2 “ % 0 2 jsin T) cos t 

trim lift coefficient f W - cos M 

rolling-moment coefficient ^ ^° — - ^ ■ D ^^H Qmen ^ j 

yawing -moment coefficient [ m Q men t ) 

\ qSb ) 

lateral-force coefficient ( Lateral i . r _ orc e \ 

V qS I 
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T l/2 


yawing -moment constant 
time, seconds 

nondimensional time parameter based on span (Vt/b) 
differential operator 

operator in Laplace transformation 

time for amplitude of oscillation to damp to one-half its 
original value, seconds 



The subscript o is used to indicate initial conditions and a bar 
is used to denote variables in the operational equations* 
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ANALYSIS 

Nonlinear Stability Derivatives 


The assumptions made with regard to the nonlinearity of the stability 
derivatives Cnp and C nr are shown in figure 1. For all three cases 

presented in the figure, C n is equal to zero for -2° < (3 < 2°, a region 
which is subsequently referred to as a dead spot. Thus when the airplane 
is within the dead spot, the value of the directional stability deriva- 
tive C n p is zero. Since the damping-in-yawing derivative C nr is a 

direct function of C n p contributed by the tail, C nr was also assumed 

to be zero for values of -2° < p < 2°. In the region outside of the 
dead spot, each one of the cases represents a different type of variation 
of C n with p in order to simulate the effect of several possible flow 
conditions on the side force acting on the vertical surface . For cases 1 
and 2, C n p = 0.28 and for case 3> C n p = 0.4l. The corresponding value 

of C nr for all three cases is -0.39- It should be noted in figure 1 

that for cases 2 and 3> C n = 0 at p of 2° and -2°, whereas for 
case 1, C n has a finite value at p of 2° and -2°. 


Method of Calculating Motion 

Since the nonlinearities shown in figure 1 can be treated as linear 
derivatives of different values within and outside of the dead spot, the 
airplane motion is calculated on the basis of classical linear theory. 
The equations of motion and the general method of calculating the motion 
of an airplane are given in references 1 and 2. The methods of refer- 
ences 1 and 2 are based on the Laplace transformation which inherently 
takes into account the initial conditions of the problem. Because the 
Laplace transformation considers the initial displacements and initial 
velocities of the problem, this method is directly applicable to the 
calculation of the motion of an airplane which has nonlinear derivatives 
similar to the derivatives presented in figure 1. 
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The nondimen sional linearized, lateral equations of motion, referred 
to the stability axes, are: 

For rolling 

2|J b ( K X 2r> b 2 0 + K XZ D b 2l|r ) = C Zp p + + £ C Z r I V !f 

for yawing 

+ ^Z°b^) - c n/ + jCnj.V + |c nr V + C nc J> (1) 

and for sideslipping 

2 ^b^ D b |3 + D b^) = c YpP + ^Yp^b^ + C L^ + ^TjPb* + ^ C L tan 

The Laplace transformation of equations ( l) , with the use of the 
symbol a for the operator, is: 

" ^Ip^? + ^b K XZ a2 " " C 2pP = 

+ ( D b0)o] - ^l/o + 2 ^b K Xzho + M°] - &l*o 

^b K XZ ff2 - ^np^ 0 + ^b K Z 2a2 " 2°^°) * ' C n p P = 

2 ^b K XZ K + Mo] - ^00 + 2 ^b K Z 2 [^o + M)°] “ 2°n r ^o + ^ 
2%° ~ C lj0 + ^b a " 1?Y t ° - C L tan ^ ^ 2^0 - Cyp^P = 

“ ipYp^o + ^b - 2°Y r ^ + 2 ^o 



M2) 


Equations (2) representjbhree simultaneous algebraic equations which can 
be solved for p, 0, ty, and their derivatives by the method of determinants. 
For example. 




£Lsl 

F(o) 


( 3 ) 
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vhere A is the characteristic lateral-stability equation 

AcA + Bct 3 + Ccr 2 + Dcr + E 


and 


= A-^cj3 + 


B 1 CT C 


+ Cqa 


+ D 


1 


The expressions for A, B, C, D, and E, in terms of the mass and aero- 
dynamic parameters of the airplane, are given on pages 27 and 28 of 
reference 1. The coefficients of the Aq equation are 


A 1 " 


Po 


8p h 3( K 


2v- 2 

•X K Z 



B i = S^o I 




^b 2c L( K X 2K Z 2 " K XZ 2 ) + ( D b^) o 2 ^b 2c Y p ( K X 2K Z 2 " K XZ 2 ) 


1 ^ |a b 2c L tan r(K x 2 K z 2 - %z 2 ) + ( D b^) o 2 ^b 2 ^ C Y r “ ^b) (% 2k: : 


2 ^b 2 (%Z C Z r - % 2c n r ) + 2 ^b 2 ( K XZ C n p “ K Z 2c i 


C n c ^‘ 7 i b K XZ C Y p + ^b K X 2c Y r " ^b^X^ 
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C 1 = S*o 


^b C L^ K XZ C Z r " K X 2c n r ^ + ^b C L^ K XZ C n p " K Z 2c Z p ) 


(M° ^ K X 2 ( C n p c Y r " C n r C Y p ) + 2 ^b % 2 (^L - ^ 


2 ^b 2 %z( C Z p ~ 2 %Z C l) + ^b K Xz( C Z r C Y p - C Z p C Y r ^ 




H b C L tan j(k xz C 1t - C^K x 2 ) + n b C L tan 7^^ - K Z 2 C^ 


( D b t )o ^b K Xz( C n p C Y r “ C n r C Y p j + ^b K Z 2 ( C Z r C Y p “ C Z p C Y r j + 


V b 2 C L tan 7 ^% 2 K z 2 - % z 2 ^ + 2 n b 2 ^K z 2 C Zp - % z C nj 


2^b ^ Z p C n r ~ C Z r C n p ^ + C nc^Yp C Z r - 2 ^b K XZ C L + 2 ^b K X 2c L tan 7 + 


^t^Zp " ]^Z p C Y r ^ 


D 1 = 


lj^L( c Zp C n r - C Z r C n T 


(^o 


M b C L tan 7(K X C n ^ - K xz C Zp '| + 


u b C L^ ] 


K XZ C Z r - *X C 


n r) 


+ t, 


ydy tan 7 /C 7 C - Cy C_ 

4 t* I Zp n r z r n p 


(V)o ^b C L( K Z 2c Z r " K XZ C n r ) + ^b C L tan ?( K XZ C n p - K Z 2 C 2] 



c\2 r 2 P 


- ^Z^L tan 7 
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The solution of equation ( 3 ) ; which will result in a time history 
of 3 as a function of s-^, is obtained from the Heaviside expansion 
theorem (reference 3): 


m 


P = 


n=l 


F'(X n ) 


^n s b 


( 4 ) 


where X n are the roots of F( cr) set equal to zero. Similar solutions 
are derived for 0, ijf, Db0, Db 1 l f , and EbP* The time scale is readily 

converted from s-^ units to t units by the equation t = ^st, • 


The values of the stability derivatives and mass characteristics 
used in the calculations are presented in table I. The table is divided 


into two columns which differ only in the values of C nQ and C : 


of 


P = 2°, 


*np ailu - °n r 

the airplane for the cases where the airplane is either outside of or within 
the dead spot. From the analytical solution of the motion, based on the 
mass and aerodynamic characteristics of the first column of table I and an 
initial condition of 0 = 5°, the time history of 3 was computed for 
several values of s-^ until the value of s-^ for which 3=2° was 
reached. For values of s^ greater than the s-b which results in 
this analytical solution is incorrect since the airplane has now entered 
into the dead spot and the values of C n p and C n ^ are zero. Thus, a 

new solution must be calculated with the use of the values given in the 
second column of table I with new initial conditions. The new initial 
conditions are determined by substituting the value of s-^ at which 
3=2° in the original analytical solutions of 0, t, Db0, Dbt, 
and D^p. Once these initial conditions are known, another set of 
analytical solutions are computed for 3, 0, 1 t f , and their derivatives 

from equations (3) and (4). This procedure is followed every time 3 
crosses through 2° or -2°. The final resultant motion in sideslip is 
the sum of all the analytical solutions in 3, each one of which is 
correct only for a particular interval of time. 


The constant Cn c is introduced into the yawing -moment equation of 

equations ( l) , since the value of the yawing-moment coefficient due to 
sideslip is C n pp + C nc for the condition of the airplane having the 

dead spot in cases 2 and 3 of figure 1. The values of C n are |o. 00977 j 
and [ 0 . Ollf-3 1 for cases 2 and 3> respectively. The sign of c Cn c is 
opposite to that of 3. For case 1 of figure 1, Cn c = 0. 
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It is apparent that the procedure employed is a time-consuming 
process and subject to the possibility of many computational errors due 
to the magnitude of the computations. The final solution can be obtained 
however, in a relatively short time through the use of automatic digital 
computing machines. 


RESULTS AND DISCUSSION 


The effect of the nonlinear stability derivatives on the lateral 
motion was investigated for the airplane described by the mass and aero- 
dynamic characteristics given in table I, with three different values for 
the damping of the lateral oscillation as calculated on the basis of 
derivatives constant with amplitude. Since the damping was varied 
arbitrarily by assuming different values for the angle of inclination of 
the principal longitudinal axis of the airplane to the flight path tj, 
three values of tj, - 2 °, 0°, and 2 °, were selected which correspond to 
a damping of the lateral oscillation, expressed in terms of ^±/ 2 , of 5 * 6 , 
3.0, and 1 . 8 , respectively. The motion of the airplane in sideslip, due 
to an initial disturbance in sideslip of 5 °; for the three values of tj 
is shown in figure 2. Since these motions are calculated on the assump- 
tion of derivatives constant with amplitude, the amplitudes of the motion 
decrease exponentially with time and will eventually reduce to zero. As 
can be noted in the first column of table I, the C n p for cases 1 and 2 

is 0.28; whereas the C n p for case 3 is O.kl. The motions presented in 
figure 2 are for C n p = 0 . 28 ; however, the motions for C n p = 0.4l would 

exhibit oscillations of approximately the same damping and a slightly 
smaller period. 

The motions of the airplane in sideslip, showing the effect of the 
nonlinearities illustrated in figure l( a) , l(b), and 1 (c), are presented 
in figures 3 to 5> respectively. In all cases, an initial disturbance 
in sideslip of 5° was assumed. The pronounced effect of the nonlinearities 
on the lateral motion is noted by a comparison of figure 2 and either 
one of figures 3* or 5* In all three figures (figs. 3 to 5) the 
motion for tj = 2 °, the most stable case, approaches a constant value. 

The analytical solution of the motion for the case of tj = 2° in figure 3 
indicates that, within the dead spot, the airplane will oscillate at a 
period of 6.56 seconds and T^yg = 3*38 seconds and will eventually 

approach the value of 0 = -0.0092°. Similar motions would be obtained 
for the case of rj = 2° in figures 4 and 5. As tj is decreased, the 
damping of the oscillatory motion depends upon the nonlinearity assumed 
and the values of 'r\. In figure 3/ the motion for tj = 0 U damps at a 
slow rate at the large amplitudes until the oscillation reaches an 
amplitude of approximately 2.^° where the damping of the oscillation is 
zero. The period of the oscillation increases from 1.5 to 1.85 seconds. 
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For the case of r\ = -2°, a very lightly damped oscillation is apparent 
within the first few seconds and the airplane may be considered to be 
neutrally stable at an amplitude of ^4,5° « In figures 4 and 5 the 
motion for tj = 0 clearly indicates that the damping is decreasing as 
the amplitude decreases and the period of the oscillation increases ; 
for t] = -2 , the oscillatory motion is slightly unstable. A neutrally 
stable oscillation would be expected to occur in figures 4 and 5 for the 
combinations of a value of T) between 0° and -2° and the dead spot 
assumed in the calculations or for r\ = -2° and a smaller dead spot* 


In general, the results indicate that the damping of the lateral 
oscillation calculated with the use of derivatives constant with amplitude 
is a determining factor in the type of motion obtained where nonlinear 
derivatives are present. As the inherent damping of the lateral 
oscillation decreases , a smaller dead spot will result in a neutrally 
stable oscillation* Obviously, if the inherent damping is zero, a 
neutrally stable oscillation already exists with zero dead spot. 


Some additional calculations were made for the case where the air- 
plane is disturbed within the dead spot. The motions for an initial 
condition of (3 = 1° were computed for r\ = -2° and 0° with the assump- 
tion of the nonlinearity described in figure l(b). The results are 
presented in figure 6. It should be noted that the only difference between 
figures 4 and 6 is the initial condition assumed in the calculations. In 
figure 6, the motion for r\ = -2° is unstable and gradually approaches 
the amplitude and period of the motion for the case of T] = -2° in 
figure 4 , The motion for rj = 0° in figure 6 is slightly unstable and 
will probably increase until its amplitude and period are in close 
agreement with the motion for the case of r\ = 0° in figure 4. Calcu- 
lations have indicated that the oscillatory motion of the airplane within 
the dead spot will double amplitude about every 4 seconds for r\ = -2° 
and about every 30 seconds for t] = 0°, If the motion' is unstable within 
the dead spot, either the airplane motion will be neutrally stable with 
an amplitude equal to or greater than the amplitude of the dead spot or 
the motion will be unstable. 


The loss in damping and the increase in period which appeared in 
some of the lateral oscillation's in figures 3 to 5 can be attributed to 
the type of nonlinearity assumed. From classical dynamic stability theory, 
it is well known that the damping of the oscillation is a function of C 


and the period of the oscillation is a function of C 


n r 


n P‘ 


If the airplane 


is considered as a mass-spring dashpot system, C n ^ is the equivalent 
spring constant of the system and C 


n corresponds to the damping 


constant contributed by the dashpot. Thus as 


Cnr 


is reduced the period 


increases and as C n is reduced the damping decreases. 
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CONCLUDING REMARKS 


The results of the investigation made to determine the effect of 
nonlinearities assumed in the analysis on the lateral stability indicate 
that under certain conditions a motion is obtained which has different 
rates of damping for the large and small amplitudes of motion, with 
very little damping at the small amplitudes. In general, the period 
of the resultant oscillation increases with time. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Va. , September 19, 1950 
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TABLE I 

STABILITY DERIVATIVES AND MASS CHARACTERISTICS 
OF THE AIRPLANE CONSIDERED IN THE ANALYSIS 


Derivative or characteristic 

Outside of dead spot 

Within dead spot 

W/S, Ib/ft 2 ..... 

80 

80 


101.1 

101.1 

P, slugs/ft3 

0.00089 

0.00089 

V, ft /sec . 

753 

753 

c L 

0.318 

0.318 

b, ft 

27.7 

27.7 


0 

0 

K Z 2 

0.0573 

0.0573 

K X 2 

O.OO69 

0.0069 

9 per radian 

P 

-0.462 

-0. 14-62 

C n p, per radian 

-0.0155 

-0.0155 

C 2^, per radian 

-0.126 

-0.126 

Cy , per radian 

0 

0 

Cy y per radian ........ 

0 

0 

d.eg . 

-2.0, 0, 2.0 

-2.0, 0, 2.0 

C n , per radian 

-0.392 

0 

C n p (cases 1 and 2), per radian 

0.28 

0 

C n p (case 3), per radian . . . 

0.4l 

0 






■2 Q 2 4 

ft; degrees 

(a) Case 1. 


-4 -a O Z 4 -4-2 O 2 4 

(3j degrees / 3 , degrees 

(b) Case 2. (c) Case 3- 


Figure 1.- Three types of nonlinear stability derivatives assumed in 

the analysis. 









t , seconds 


Figure 5.- The effect of the nonlinear derivatives described in 
figure 1(c) on the motion of an airplane. 
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t , seconds 

Figure 6.- The effect of the nonlinear derivatives described in 
figure l(b) on the motion of an airplane. Initial disturbance 
in sideslip of 1°. 
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